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In this article, I will report a Lax pair structure, a Backlund-Darboux transformation, 
and the investigation of homoclinic structures for 2D Euler equations of incompressible 
inviscid fluids. 
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1 Introduction 

The governing equation of turbulence, that we are interested in, is the incompressible 2D Navier- 
Stokes equation under periodic boundary conditions. We are particularly interested in inves- 
tigating the dynamics of 2D Navier-Stokes equation in the infinite Reynolds number limit and 
of 2D Euler equation. Our approach is different from many other studies on 2D Navier-Stokes 
equation in which one starts with Stokes equation to prove results on 2D Navier-Stokes equation 
for small Reynolds number. In our studies, we start with 2D Euler equation and view 2D Navier- 
Stokes equation for large Reynolds number as a (singular) perturbation of 2D Euler equation. 
2D Euler equation is a Hamiltonian system with infinitely many Casimirs. To understand the 
nature of turbulence, we start with investigating the hyperbolic structure of 2D Euler equation. 
We are especially interested in investigating the possible homoclinic structures. 
| In |l]], we studied a linearized 2D Euler equation at a fixed point. The linear system decou- 

ples into infinitely many one-dimensional invariant subsystems. The essential spectrum of each 
invariant subsystem is a band of continuous spectrum on the imaginary axis. Only finitely many 
of these invariant subsystems have point spectra. The point spectra can be computed through 
continued fractions. Examples show that there are indeed eigenvalues with positive and negative 
real parts. Thus, there is linear hyperbolicity. 

In |g] and Q , a Lax pair and a Backlund-Darboux transformation were found for the 2D Euler 
equation. Typically, Backlund-Darboux transformation can be used to generate homoclinic 
orbits §]. 

The 2D Euler equation can be written in the vorticity form, 

d t n + {*,n} = o , (i) 

where the bracket { , } is defined as 

{f,g} = (d x f)(d y g)-(d y f)(d x g) , 

where \E' is the stream function given by, 

u = , v = d x ^ , 

u and v are the velocity components, and the relation between vorticity £1 and stream function 
* is, 

fl = d x v - d v u = . 
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2 A Lax Pair and a Darboux Transformation 

Theorem 1 (Li, [§]). The Lax pair of the 2D Euler equation ^) is given as 

L(p = \ip , , , 

d t <p + A<p = , y ' 

where 

Lip = {n,ip} , Aip = {y,<p} , 

and X is a complex constant, and <p is a complex-valued function. 
Consider the Lax pair @ at A = 0, i.e. 

{<n, P } = o, (3) 

dtP + {V,p} = 0, (4) 
where we replaced the notation ip by p. 

Theorem 2 (Li and Yurov, Q). Let f = f(t,x,y) be any fixed solution to the system (||, 
we define the Gauge transform Gf: 

P = GfP = 75- \Px ~ (d x In f)p] , (5) 

and the transforms of the potentials and Vl/ : 

^> = ^> + F,n = n + AF, (6) 
where F is subject to the constraints 

{U,AF} = 0, {n + AF,F} = 0. (7) 

Then p solves the system (B, at (Cl,^). Thus and (^) form the Darboux transformation 
for the 2D Euler equation (yj and its Lax pair (0, 

3 Preliminaries on Linearized 2D Euler Equation 

We consider the two-dimensional incompressible Euler equation written in vorticity form([l]) 
under periodic boundary conditions in both x and y directions with period 2ir. We also require 
that both u and v have means zero, 

27T f2n e-l'K r2n 

u dxdy = v dxdy = 0. 



>o J Jo Jo 

We expand £1 into Fourier series, 

fi= ^e ihX , 
kez 2 /{o} 

where = uJ^ , k = (fci,^) 7 ', X = {x,y) T . In this paper, we confuse with (0, 0) T , the 
context will always make it clear. By the relation between vorticity Q and stream function vl/, 
the system (|l]) can be rewritten as the following kinetic system, 

UJ k = ^ A (P> "?) U V^q , (8) 
k=p+q 
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where A(p, q) is given by, 

Mp,q) = \[\q\~ 2 - W\piq2 -P2Qi) , (9) 

where \q\ 2 = q\ + q\ for q = ((?i, (?2) T > similarly for p. 

We denote {^k}kez 2 /{o} by u. For any fixed p £ Z 2 /{0}, we consider the simple fixed point 

uj*: 

u* p = T, u* k = 0, if k + p or -p, (10) 

of the 2D Euler equation (||), where T is an arbitrary complex constant. The linearized two- 
dimensional Euler equation at u>* is given by, 

Wfc = A(p, k - p) T uj k _ p + A(-p, k + p)T uj k+p . (11) 
Definition 1 (Classes). For any k G Z 2 /{0}, we define the class to be the subset of Z 2 /{0}: 




Figure 1. An illustration of the classes E^, and the disk D\ p \- 



See Figure p] for an illustration. According to the classification defined in Definition [D, the 
linearized two-dimensional Euler equation ([ll]) decouples into infinitely many invariant subsys- 
tems: 

"k+np = A (P: k + {n- l)p) T U} k+(n _ 1)p 

(12) 

+ A(-p, k + (n + l)p) f u %Hn+1)p . 
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Theorem 3. The eigenvalues of the linear operator Ct defined by the right hand side of (f?H), 
are of four types: real pairs (c,—c), purely imaginary pairs (id,— id), quadruples (±c±id), and 
zero eigenvalues. 

The eigenvalues can be computed through continued fractions. 

Definition 2 (The Disk). The disk of radius \p\ in Z 2 / {0}, denoted by D\ p \, is defined as 



D ]p] = {k£Z 2 /{0} 



\k\<\p\ • 



Theorem 4 (The Spectral Theorem). We have the following claims on the spectra of the 
linear operator : 

1. If Yi^C\ D\ p \ = ; then the entire £2 spectrum of the linear operator C^, is its continuous 



spectrum. See Figure^ where b = — i|r||p|" 



Pi h 
P2 h 



2. If X^, n -D| p | 7^ ; then the entire essential £2 spectrum of the linear operator Cf, is its 
continuous spectrum. That is, the residual spectrum of is empty, a r (Cf,) = 0. The 
point spectrum of is symmetric with respect to both real and imaginary axes. See 
Figure ||. 
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Figure 2. The spectrum of Cf,. 
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4 A Galerkin Truncation 

To simplify our study, we study only the case when u>k is real, VA: G Z 2 /{0}, i.e. we only study 
the cosine transform of the vorticity, 

= oik cos(k ■ X) , 

fcez 2 /{o} 

and the 2D Euler equation (|l|||) preserves the cosine transform. To further simplify our study, 
we will study a concrete line of fixed points (^) with the mode p = (1, 1) T parametrized by T. 
When r / 0, each fixed point has 4 eigenvalues which form a quadruple. These four eigenvalues 
appear in the only unstable invariant linear subsystem labeled by k = (—3, — 2) T . See Figure 
H for an illustration. We computed the eigenvalues through continued fractions, one of them is 
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Figure 3. The collocation of the modes in the Galerkin truncation. 



A = 2A/|r| = 0.24822302478255 + i 0.35172076526520 . (13) 

We hope that a Galerkin truncation with a small number of modes including those inside the 
disk D\ p \ can capture the eigenvalues. We propose the Galerkin truncation to the linear system 
flj~2| ) with the four modes k + p, k + 2p, k + 3p, and k + 4p, 

UJ 3 = A 2 TU!2 — A^Fui , 

Co 4 = AzTtu-s . 
From now on, the abbreviated notations, 



A n = A(p, k + np), A n 



A(k + mp, k + np) , 



(14) 
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will be used. The eigenvalues of this four dimensional system can be easily calculated. It turns 
out that this system has a quadruple of eigenvalues: 



A = ±-^=Jl±iV35= ± ( M x 0.7746 x e ±idl , (15) 



2VW v V 2 

where 6± = arctan(0.845), in comparison with the quadruple of eigenvalues (fb|), where 



A=± 



Q x 0.43 x e ±id2 , 



and 62 = arctan(1.418). Thus, the quadruple of eigenvalues of the original system is recovered 
by the four-mode truncation. We further study the corresponding Galerkin truncation of 2D 
Euler equation: 

Cj\ = —A 2 U! p L02 , 

L02 = A\ UJ p U)\ - A 2 UJ p UJ 3 , 

uj 3 = A 2 ujp uj2 - Ai ujp uii , (16) 

w 4 = A 2 ujp uj 3 , 

UJp = Ai >2 (U> 3 W 4 - UJi U> 2 ) , 

and the equations for the decoupled variables ujq and ui§ are given by, 

= ~A\ ujp ui , 
ti;5 = A\ ujp UJ4. . 

where 

3 1 

Ai = -— , A 2 = - , A 3 = A 2 , A 4 = Ax , 

4 

Ai >2 =A 1 -A 2 = —=, A 2 , 3 = , A 3A = -A lt2 ; 
5 

There are three invariants for the system (|l~6|): 

I = 2Ai^(uj 1 uj3 + ^2^^) + A 2 ujp , (17) 
U = A^uol + uj\) + A 2 {uj 2 2 + ujI) , (18) 
J = u? v + U)\ + U?2 + w| + ujj . (19) 

J is the enstrophy, and U is a linear combination of the kinetic energy and the enstrophy. / is 
an extra invariant which is peculiar to this invariant subsystem. With I, the explicit formula 
for the hyperbolic structure can be computed. 

The common level set of these three invariants which is connected to the fixed point (|i"o| ) 
determines the stable and unstable manifolds of the fixed point and its negative — uj*: 

ujp = -T,uJ n = (neZ). (20) 

Using the polar coordinates: 



uj\ = r cos 9 , UJ4 = r sin 8 ; LO2 = p cos $ , uj 3 = p sin 1? ; 
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we have the following explicit expressions for the stable and unstable manifolds of the fixed point 
dlOD and its negative (p0|) represented through the homoclinic orbits asymptotic to the line of 
fixed points: 



r tanh r . 



A 2 -A l 
A 2 



T sech t , 



2k 



In cosh r + 6q , 



(21) 



A 2 



+ 1? 



arcsm 



1 / ^2 

2 A/ -Ai 



7r + arcsin 



1 ./ M 

2 V -Ai 



(k > 0) , 
(k < 0) , 



where Ax and A 2 are given in (|i6[), r = nYt + ro, (to, 0q) are the two parameters parametrizing 
the two-dimensional stable (unstable) manifold, and 



k = \j-A x A 2 cos{9 + &) = ± V / -A 1 A 2 \ 1 + 



Ao 



4Ai ' 

The two auxilliary variables loq and oj^ have the expressions: 



LO 



where 



a 



aft 
1 + /? 2 



-AiTkT 1 



sech r sin[/3 In cosh r + 8q] — - cos[/3 In cosh t + 9q\ 

sech r •{ cos[/3 In cosh r + #o] + -5 sin[/3 In cosh r + 9q \ 

P 



An 



2k 



A 2 -A 1 

The graphs of these homoclinic orbits are spirals on a 2D ellipsoid, with turning points. 



5 Conclusion 

Certain newly developed results on 2D Euler equation have been discussed, which include a Lax 
pair, a Darboux transformation, and the investigation on homoclinic structures. 
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